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Presently used substrates for integrated circuit chip mounting and 
interconnection provide limited heat sink capability and preclude the 
effective use of matched impedance transmission line interconnection. 
These characteristics can be improved by the use of very narrow microstrip 
lines on thin dielectric layers backed by thermally and electrically con- 
ductive material. We propose a model for calculating the effective thermal 
resistance of layered substrates and make calculations for a variety of 
materials. The thermal problem is considered for heat extraction either 
through the edge or through the face of the substrate. Improvements in sub- 
strate thermal conductance of an order of magnitude appear realizable 
compared to presently used substrates (for example, 0.625 mm alumina). 
We analyze the electrical parameters of narrow (0.02-0.1 mm) microstrip 
lines on coated substrates in the frequency range 10 MHz to 10 GHz. The 
characteristic impedance, propagation delay time, and attenuation are 
found to be frequency dependent, and efforts to minimize this frequency 
dependence by magnetic loading result in greater, though more constant, 
delay times. Losses may be significant (~2-5 dB/cm in many cases), 
but short line lengths due to dense circuit packing allowed by improved 
heat dissipating capability will minimize this disadvantage. 

I. INTRODUCTION 

A necessary component in realizing the potential of integrated circuit 
technology is the substrate upon which the integrated circuit chip is 
mounted. The substrate must meet several requirements to avoid 
deterioration of the system's performance: 

(i) The substrate must act as a rigid, reliable mounting platform for 
the chip. 
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(ii) The substrate must allow the dissipation of heat produced on 
the chip without causing an increase in temperature which would detri- 
mentally effect circuit operation. 

(Hi) The substrate must permit electrical connections to and inter- 
connections between chips in a manner compatible with desired circuit 
performance. 

(iv) The substrate fabrication and chip mounting must be realizable 
at a cost consistent with the cost requirements of the system. 

Presently used substrates (for example, 0.625 mm thick A1 2 3 ceramic 
wafers) are limited in their attainment of objectives ii and Hi 
outlined above. Heat sink capability is limited by the low thermal con- 
ductivity and relatively large thickness of material necessary for mechan- 
ical rigidity and by the difficulty in extracting heat from the substrate 
itself. High frequency (>500 MHz) performance of beam-leaded in- 
tegrated circuits suffers from the transmission line mismatch necessary 
with conventional substrates as is discussed below. 

A new materials configuration, as shown in Fig. 1, is envisioned which 
would simultaneously satisfy the four basic substrate requirements 
listed above. A thin dielectric layer (0.001 mm to 0.05 mm thick) 
deposited or grown on a massive electrically and thermally conductive 
ground plane could allow efficient heat dissipation, effective electrical 
interconnection and provide a rigid, inexpensive bonding plane. It is 
the purpose of this paper to calculate the thermal and electrical prop- 
erties of the coated substrate configuration suggested here. 
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Fig. 1 — Schematic drawing of layered substrate, 
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II. THERMAL PROPERTIES 

The model used to represent the proposed substrate for calculation of 
thermal properties is shown in Fig. 2. All constants may be arbitrarily 
set, and the heat flow may be subjected to either of two distinct bound- 
ary conditions: the sides of the disk may be held fixed at a constant 
temperature and the top and bottom faces insulated, or the bottom face 
may be held at a constant temperature and the sides and top insulated. 
These cases closely approximate two heat sinking configurations fre- 
quently found in practice. The cylindrical geometry was chosen for 
ease of calculation (the square or rectangular problem, although alge- 
braically simpler, results in a double series which presents formidable 
convergence problems when one is interested in a numerical result). The 
specific geometry should have little effect on the calculated values, and 
any reasonable geometry (held constant throughout the calculations) 
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Fig. 2 — Model used in calculations of thermal properties. 
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should provide a reliable comparison between alternate substrates. 
The heat input was assumed uniform over a circle of radius A, although 
other heat input distributions might have been chosen. As in the choice 
of geometry, this should have a minimal effect on the results. (Currently 
used beam leaded chips distribute heat over an area larger than the 
silicon dimensions. It can be assumed that the entire chip-beam lead 
structure thus forms a relatively uniform heat source on a substrate. 1 ) 

2.1 Algebraic Solution 

The thermal problem has already been solved for a homogeneous 
cylinder. 2 In the case of a layered cylinder, it is necessary to solve La- 
place's equation simultaneously in both the upper and lower regions 
and match boundary conditions at the interface. In the calculations 
described in this paper, no discrete thermal resistance was inserted 
at the interface so the matching of boundary conditions reduced to 

where K x , T x , K 2 , T 2 are the thermal conductivities and temperatures 
in regions 1 (upper) and 2 (lower). The insertion of a thermal resistance 
at the interface is trivial but was approximated as zero due to the inti- 
mate interface contact obtained by the method of materials preparation 
envisioned for these substrates. 

Laplace's equation was solved for the two sets of boundary conditions 
(sides grounded and face grounded) by the method of separation of 
variables. In cylindrical coordinates, V 2 T = reduces to 

T = * + *i2 

+ EI * m . n [A m J n (8 n r) + 5 m AT m (5 n r)]{ S n m md\ 

[exp (dsM)] (2) 

where boundary conditions exist to impose eigenfunction solutions in the 
r direction, and J m and N m are Bessel functions of the first and second 
kind respectively. The boundary condition on the grounded surfaces was 
given as T = while dT/d (normal) was set equal to zero on the insulated 
surfaces. On the top surface (z = 0) the boundary condition was written: 



\=l 



dT, ) T L 0<r<A, 2 = 0; ._. 

~dz ~ ' (3) 

1 A <r <B, 2 = 0. 
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/ is the value of heat input in (watts/unit area). This boundary con- 
dition was met by using the orthogonality properties of the Bessel 
functions in a Fourier-Bessel series as outlined in Morse and Feshbach 
and other texts. 3 

For the case in which the sides are grounded at T = 0, the tempera- 
ture distribution is found to be 



i M*J|Mgi 



c, cosh hg 



£-)+ C2 smh(| 2 



where /3„ is the nth zero of J Q (r) and 

m tanh qj : te : v tanh v-rJ tanh l g J - 



(4) 



- (*?)[' " 1 " - Pfi - FiF)] + I - FlF) ' 



for <z < W> ; 



ta„h^-fM-: 



(5) 



tanh (-^1 - ^ - tanh [-^j tanh (-^ ijj + ^tanh ^ ' g "J 

for IF, < z < TF, + »T, ; 
1 , for < z < W, ; 



tanh ^__j^ _ _ _ tanh ( ^-j tanh [-^ *]J + ^ tanh ^ ' g " J 



,(6) 



for ir, < z < if, + if, . 



For the case in which the bottom face is grounded at T = 0, the 
temperature distribution is found to be 

T = 4 /(di + d 2 Z) 



+ z 



_ i ^^-|V«(«- r> 



where a„ is the nth zero of Jj(r) and 



g(a.) [* C ° Sh (f ) + ^ Shlh (¥. 



d, - 



K, 



< z < TF, 



TFx < z < W\ + IF, 



(7) 



(8) 
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< Z < W t 

} 

FT, < z < W x + W 2 



& = . 



if, 



(9) 



rf 3 = ■ 



for < z < W, ; 

^(^±B)[l-^(f)] 



(10) 



rf. = 






for W, < z < II', + ir 5 ; 



fl for < 2 < 1)', 



tanh 



tanh (^-g-J + ^ - lj tanh ^-g-J tanh [ g J - ^ 

for IF, < z < W, + FT, 



(11 



A quantitative measure of the power dissipating capabilities of a 
substrate is given by the thermal resistance 



„ /de grees C 

# the r in I 77 

\ watt 



AT(°C) 



T 

-* max 



q (watts input) irA j 
T max clearly occurs at z = 0, r = so 



(12) 



-Rihrm. (side grounded) = 



££ 



,44) 



K,A ^pl J»(pJ 
_ tanh ^__j - ^ - lj tanh \ rr ) tanh ^-^— - 1) - 1 



-»)t-s--(4 k )-( Mk ^]+g-T s # s 5 (131 



and 



JC, , 2 ~] • / '( a "B) 



W,K, + IJVv, ^ v * 



Biii.™ (face grounded) irB 2 K x K ; T rf,.'l r<>' •/;(«.) 



(14) 
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2.2 Computer Program 

A Fortran IV computer program was written for an IBM 360-50 to 
evaluate i2th.rm-.id. and R tberrn ., ace • The first forty values of alpha 
and beta (the zeros of Jj and J ) are inserted in the program from the 
tabulated values of Watson. 4 The values of the higher zeros are obtained 
in the program by analysis of the asymptotic approximations for J t 

and Jo • 

Since the expressions for R thl:Tm take the form of infinite series, no 
completely precise evaluation can be made. Figure 3 shows the ap- 
proximation to ftthcrm-fao. given by truncating the summation after 
varying numbers of terms. Plotting data similar to Figure 3 for a range 
of all input parameters (A, B, K X ,K 2 ,W X , W 2 ) disclosed several trends 
in the convergence of the series. The "period" of the oscillation was 
inversely correlated to the ratio A/B. This is mathematically consistent 
with the use of a Fourier-Bessel series to describe the top face of the 
disk. For ratios of A/B less than 0.02, the first maxima in the expression 
for # th erm-.id. is not reached until partial sums with greater than 300 
terms are evaluated. The amplitude of the oscillations in the function 
ftth.rm (number of terms) was positively correlated to K 2 /K x and Wi/W 2 ■ 
The damping of the oscillation (per cycle) appeared to be relatively 



2.4 








P.O 




o 


<P°o 


'OOOCX*^ 


X 
111 

Q. 


o 






.°1.2 

z 


o 






in 
A?" 0.8 


o 
o 






0.4 


o 










o 

> 


1 1 1 1 


1 1 



100 150 200 250 300 

NUMBER OF TERMS IN SUMMATION 



Fig. 3- Partial sums approximating flf ace for A = 0.0625, B = 3.75, K t = 0.20, 
A' 2 = 2.05, W: = 0.0005, W« = 0.0625. 
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independent of input data. These variabilities in the behavior of .Rtherm 
made determination of a precise value difficult. Additionally, only the 
first 314 terms of the series could be evaluated. The operation of the 
program was terminated at this point by a floating point overflow which 
occurred in the IBM-SSP subroutine used to generate the Bessel func- 
tions. Since the largest arguments of the Bessel functions are intrinsic 
to the problem, and are not functions of the input parameters, nothing 
could be done to circumvent this difficulty. It thus became necessary 
to evaluate R t h er m on the basis of a truncated series of 314 terms of less. 
A graphical outline of the method used to extract an approximation to 
the asymptotic value of i2 therm from a finite number of terms in shown 
in Fig. 4. The first two local maxima, M x and M 3 , and the first local 
minimum, M 2 , were recorded; and the averages, Ai and A. t , were cal- 
culated where 



A, = 



A, = 



Jfi 


+ M 2 


M 2 
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(15) 
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Fig. 4 — Pictorial of extremum averaging technique used for improving series 
convergence. 
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In Appendix A, a theorem is stated and proved to show that the asymp- 
totic value of #th„rm is less than A t and greater than A 2 . An average of 
the two, (A x + A 2 )/2, was used as a most probable value of R theim . 
As a test of the accuracy of the derivation of R theTma - a id<, and fl t herm-faoe 
and the computer programs for generating these values, K 2 was set 
equal to K y (equivalent to a uniform, single layer cylinder) and values 
of #tb.r„, computed for W x = 0.625 mm, W 2 = 0; W, = 0.375 mm, 
W 2 = 0.25 mm; and W t = 0, W a = 0.625 mm. The three cases were 
in complete mutual agreement, and agreed to the accuracy given with 
the published results of D. P. Kennedy for heat conduction in a homo- 
geneous, isotropic cylinder. 

2.3 Numerical Evaluation 

Several materials combinations were chosen for evaluation as possible 
substrates. Some of the considerations for inclusion were compatibility 
with integrated circuit techniques (for example, silicon), smooth surfaces 
for microdeposition of conductor paths (for example, polished silicon, 
cold-rolled aluminum), low dielectric constant (for example, SiON), 
and thermal expansion coefficient. The materials considered and the 
associated values of K used with them are shown in Table I. A = 0.625 
mm and B = 3.75 cm were chosen as representative of typical inte- 
grated circuit chips and substrate sizes and were used in most of the 
calculations. 

The results of these calculations for a variety of layered substrates 
are shown in Figs. 5 through 10. With high conductivity dielectrics 
(for example, A1 2 3 ) the thermal resistivity of the ground plane con- 
tributes significantly to the total thermal resistance, but resistances as 
low as that of beryllia are attainable for a considerable range of dielectric 
thicknesses. Low conductivity, glassy (for example, SiON, Si0 2 ) dielec- 
trics, can result in thermal resistances equal to or greater than that of 
alumina if glassy layers of greater than 0.03 mm thick are used. 

Table I — Thermal Conductivity of Substrate Materials 

°C) 



Material 


K (watts /cm 


Al 


2.05 


A1 2 3 (deposited) 


0.20 


AUO3 (sintered) 
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0.012 


SiON 


0.012 
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Fig. 5 — Aluminum oxide on 0.25 mm aluminum. 



For comparison, thermal resistances of presently available homogene- 
ous substrates were calculated and are given in Table II. 

Figure 11 shows the effect of substrate radius on /2 t u Br m- B id e (the 
effect on R t h Br m-f*c e was <2% for B ^ 0.5 cm). In the range shown, 
heat is presumably flowing radially at r = B so increases in B add to 
^tberm at a rate proportional to 1/B. Consequently, a "critical radius" 
exists above which fl t herm-side is relatively constant. 

Radiative and convective heat losses will be negligible for the high 
conductivity layered substrates discussed in this paper. A lumped value 
of convective and radiation coefficient of 0.003 watts/cm. 2 — °C 
results in an equivalent parallel thermal resistance of greater than 50°C/ 
watt for all cases considered here. 

2.4 Thermal Results 

Layered substrates consisting of dielectric layers backed by thermally 
conducting ground planes can provide substantial improvements in heat 
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sink capability over alumina ceramic substrates. Ceramic substrates 
are commonly used in 0.625 mm thickness to provide needed strengths. 
Metal ground planes could be reduced to 0.25 mm thickness while 
maintaining physical strength greater than 0.625 mm ceramic. For 
situations where heat sinking the back of the substrate is possible, 
A1 2 3 on 0.25 mm aluminum can provide up to twice the heat sink 
capability of 0.625 mm beryllia. Metal substrates have the additional 
advantage of convenient, low thermal barrier heat sink mounting or 
even direct incorporation into the heat sink structure (for example, 
heat pipes, finned substrates). 

Although the above values of fi lhorm are significant, as given, for 
relative comparison of substrate materials and geometries, their ab- 
solute significance can only be assessed in relation to other parameters 
intrinsic to an integrated circuit heat dissipating system. V. E. Holt 
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Fig. 6 — Aluminum oxide on 0.625 mm aluminum, 
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Fig. 7 — Silicon oxynitride on 0.25 mm aluminum. 

has found thermal barriers of 32°C/watt for 1.25 mm chips conducting 
heat to a substrate by means of beam leads and a resin bonding layer. 1 
This value is reduced to 3 to 5°C/watt by AuSi eutectic bonding. 
Imperfections in the heat sink will have to be assessed in each situation 
encountered. 



III. ELECTRICAL PROPERTIES 

3.1 Microstrip Transmission Lines fc 

Electrical interconnections on the proposed substrates would best be 
made in the form of standard microstrip transmission lines. 5 They can 
provide transmission line interconnects necessary for high speed pulse 
and microwave electronics while maintaining low values of crosstalk 
between adjacent lines and still provide dc and low frequency paths 
for other circuit requirements. The primary incompatibility of presently 
used substrates and beam-leaded integrated circuit chips is due to the 
impedance mismatch between the narrow (0.02 to 0.05 mm) beam leads 
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and the physically large transmission line (typically 5012 implies 0.625 
mm line width on 0.625 mm alumina). 5 The removal of this discontinuity 
would require narrow conductors deposited on a dielectric layer of 
appropriate thickness to maintain the desired transmission line im- 
pedance. The physical configuration contemplated for small microstrip 
lines is outlined in Fig. 12. To minimize the impedance discontinuity 
at the interconnection with beam leaded chips, W is restricted to 0.1 
mm maximum. To minimize crosstalk, this would suggest t d ^ 0.05 mm'. 
These figures, coupled with the electrical parameters of dielectric materi- 
als sufficiently smooth to allow the reliable deposition of 0.02 mm 
conducting paths, are consistent with the desire to develop lines with a 
characteristic impedance of 5012, the most commonly used interface 
impedance for circuitry operating above 500 MHz. Under certain con- 
ditions, as seen below, it is desirable to raise the series inductance of the 
transmission line. For this reason, provision has been made for placing a 
thin layer of high permeability material between the conducting strip 
and the dielectric. In practice the adhesion metal of a sandwiched inter- 
connect (dielectric— adhesion metal— gold) might serve this purpose. 
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Fig. 8 — Silicon oxynitride on 0.625 mm aluminum. 
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3.2 Transmission Line Parameters 

Although many transmission modes are necessarily present in an 
unsymmetrical transmission line such as an unshielded microstrip, the 
TEM mode will be considered of dominant interest in the frequency 
range considered here and the line parameters will be calculated for 
this mode. Several authors (for example, Assadourian and Rimai, 7 and 
H. A. Wheeler 8 ) have derived expressions for the characteristic imped- 
ance and other pertinent parameters for microstrip transmission lines. 
Although these formulae have been generally verified in their range of 
approximation 5,0,10 the physical scale of lines contemplated in the last 
section requires the inclusion of several terms which had been ignored 
in the previous derivations. 

The characteristic impedance, Z , and propagation constant, 7, of 
any power transmission system are given by 11 



„ _ ( R + joL V 
*° = \G + juC/ ' 

y = [(R + jcoLXG + j<oC)] J . 



(17) 
(18) 
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Ignoring losses (R, G 
written 



0), the propagation delay time can then be 



r = — ** = (LC) h . 



(19) 



The inductance of the line is given by two components: the inductance 
due to magnetic field energy storage in the strip conductor and the 
ground plane (the internal inductance); and the inductance due to mag- 
netic field energy storage between the strip conductor and the ground 
plane (the external inductance). The internal inductance of the strip 
conductor can be derived from Ramo and Whinnery 12 and is given by 
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Fig. 10 — Silicon dioxide on 0.625 mm silicon. 
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Table II — Thermal Resistance of Conventional Substrates 





'Ithorm-sido 


R tharm- face 




(degrees C/watt) 


Alumina 






0.25 mm 


103 


6.8 


0.625 mm 


46 


12 


Beryllia 






0.25 mm 


15.3 


1.02 


0.625 mm 


6.8 


1.8 



skin depth of the strip conductor at the frequency, «, of interest. [5 = 
(2p/w/*)*]. For the frequency range considered in this paper (/ > 10 
MHz) current spreading in the ground plane can be ignored for lines as 
narrow as 0.02 mm if the ground material has a resistivity ^10~ 6 
ohm-cm. For higher resistivity materials this approximation is less 
valid, especially at the lower frequencies, but the approximation will 
be made for convenience. A more rigorous analysis in the higher resis- 
tivity case would give increased internal inductance and increased 
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Fig. 11 — Thermal resistance with side grounded A = 0.625 mm Wi = 0.0125 
mm, W2 = 0.625 mm. 
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■pig. 12 — Schematic drawing of microstrip transmission lines on layered substrate. 

crosstalk through magnetic fields generated by spreading currents in 
the ground plane. The internal inductance per unit length of line due 
to fields in the ground plane can then be given by a formula similar to 
(20) 



U = 



Wu 



p. 8iph (^)- Bin (f 
«- (2t, 



cosh 



— cos 



2/, 



«„/- 



(21) 



where p, and 8„ are the resistivity and skin depth of the ground plane. 
For the dimensional range of interest (W > 2t d ) the external inductance 
per unit length is approximately given by English and McNichol 13 as 



L Mt - £ ft + HaQ 



(22) 



where n R is the relative permeability of the magnetic material. Total 
L is given by the sum 

L = L. + L„ + L ext . (23) 

The total capacitance per unit length between the strip conductor and 
the ground plane is the sum of four components 

C = C T + C + C, + C t . (24) 
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C T is the capacitance per unit length between the top of the strip and 
the ground at infinity 14 . 



C T = — In 2. 



(25) 



C„ is the capacitance per unit length, neglecting fringing, between the 
face of the strip and the ground plane. 



(-'a ~ Co e fi 



W 



(26) 



where e R is the relative dielectric constant of the insulating layer. The 
additional capacitance per unit length due to the fringing fields at the 
edges of the strip is given by Joines 15 

1 + II — exp 



C, = ^ In 



2L 



(27) 



For the case where the strip conductor is not of negligible thickness, an 
additional contribution to the capacitance is due the capacitance be- 
tween the edges of the strip and the ground plane. The expression for 
this term, derived in Appendix B, is 

J2t d + 2/, + [2(t d + /,)(/„ + 2t.)f 
L 



C. = ^ In 



(28) 



Losses on these microstrip lines will be considered as a perturbation 
in the calculation of the other parameters although this approximation is 
quite poor in several of the cases considered below. Losses in the dielectric 
will be ignored, being much less than the conductor losses for the geo- 
metries discussed in this paper. Under considerations similar to those 
for deriving the internal inductance, the equivalent series resistance is 
given by 11,12 



R = W 



smh(^)+sm [f 



^sinh^+sin^ 



_ 8 ' cosh (f) - cos (f) )B cosh (f) - COS (|*)_ 



(29) 



The attenuation per unit length of line is given by" 

a = 9^"( ne P erS )- 



(30) 



Although this expression is rigorously valid only for | R/u \ « | L |, it 
will be used in all computations made in this paper. In the case where 
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| R/co | > | L |, the line will have a significant reactive component in its 
characteristic impedance. Equation (30) for a will then contain a con- 
tribution for phase shift as well as attenuation. In this paper \Z \, 
Re (Z n ), and Re (a) will be calculated. 

The equations for crosstalk between adjacent, parallel strip lines 
derived by Kordos apply directly to microstrip lines of the dimensions 
considered here if the characteristic impedance is calculated by the 
means outlined above. In these cases, Kordos derives an expression for 
the near end crosstalk of properly terminated parallel lines of length, I. 

30 In [1 + (tJWff f1 . m /Q1 x 

Re(Z,)(«y [l-exp(-2srZ)]. (31) 

V n ,(s) and V„(s) are the Laplace transforms of the near-end crosstalk 
voltage and the input signal respectively. The far-end crosstalk is 
approximately zero. This formula is only valid for t m = 0. For t m > 0, 
V ne would be further reduced by the "keeper" effect of the magnetic 
material on the magnetic field of the driven line. 

3.3 Numerical Evaluation 

The formulas given above have been applied to layered substrates of 
the type discussed earlier. The calculations assumed that the strip con- 
ductors (Au — 2.5 X 10~°&— cm.) should exhibit a resistance of 
<0.004fl/sq. to give a dc resistance <1.5S2/cm for the narrowest line 
considered, 0.02 mm. This results in t, = 0.00625 mm. Silicon and 
aluminum were considered as possible ground plane materials. The 
maximum doping level in silicon which allows the preservation of a 
smooth surface after processing is approximately 2 X 10 1 atoms/cm 
which results in a bulk resistivity of approximately 2.5 X 10~ 3 fi— cm. 
The constant parameters used in the calculations are shown in Table III. 

The expressions given in Section 3.2 were evaluated for a spectrum 
of frequencies between 10 MHz and 10 GHz for several physically 
realizable dielectric-ground plane materials combinations with varying 

Table III — Electrical Parameters for Stripline 
Calculations 









(All dimensions 


in cm) 




p. 

Pa 

Pa 

Mr 

t. 


(Au) 
(Al) 
(Si) 


= 2.5 
= 2.5 
= 2.5 
= 100 
= 6.2. 


x io-« 

x io-° 

X 10~ 3 

> x io-" 




= 6.25 X 10" 2 
(sintered A1 2 3 ) = 
(deposited AI5O3) 
(SiON) =3.8 
(SiO : ) =3.6 


9.6 
= 9.0 
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thicknesses of magnetic material. A representative sampling of the 
results is shown in Figs. 13 through 20. These data are not meant to be 
comprehensive or complete but merely to represent the effects and 
trends observed in the calculations. 

Figures 13 through 16 give Re (Z ) and t as a function of dielectric 
thickness and/or line width. The effect of the large internal inductance of 
Si on Z and t is especially apparent for thin dielectrics (where the ex- 
ternal inductance is small). It is noted in Figs. 13 and 14 that r passes 
through a local minimum at t d ~ 0.01 mm mils. Below this value (for 
aluminum ground planes) the constant internal inductance is greater 
than the external inductance; and the capacitance, which decreases 
with t d , causes r to decrease with t d . Above the minimum, as the in- 
creasing external inductance predominates, the fringing capacitance 
does not decrease linearly with t d so t rises. As t d increases further, the 
equations used in this paper cease to be valid and the inductance will 
not increase linearly with t d . 

In Figs. 17 and 18, the solid lines show the frequency dependence of 
the characteristic impedance. The variation of Re (Z ) with oj is particu- 
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0.20 



0.16 



0.08 



0.04 



0.005 0.010 0.015 0.020 0.025 

DIELECTRIC THICKNESS IN mm 



0.050 



Fig. 13 — Impedance and delay time at 1 GHz for aluminum oxide dielectric on 
aluminum ground plane. 
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50 




0.24 



0.20 



0.16 



0.12 



0.08 



0.005 0.010 0.015 0.020 0.025 

TOTAL DIELECTRIC THICKNESS IN mm 



0.030 



Fi^. 14 — Impedance and delay time at 1 GHz for Bilicon oxynitride dielectric on 
aluminum ground plane. 



larly large in the case of Si since the frequency independent external 
inductance is small compared to the internal inductance. For pulse 
applications in which an impedance match is desirable over a broad 
frequency range, the situation may be improved somewhat by increasing 
the frequency independent component of the inductance. This may be 
done by loading the line with magnetic material, the results of which 
are shown by the dashed lines in Figs. 17 and 18. The open triangles 
show the unavoidably greater propagation delay times associated with 
the loaded lines. It can be seen from the deriving formulas that only 
the product /x R t m was considered significant, so materials of other 
permeabilities and appropriate thicknesses would be equally effective. 
The losses associated with these lines have been calculated, and a sam- 
ple of the results is shown in Fig. 19. Attenuation (dB/cm) is plotted for 
two materials configurations at each of two frequencies for a range of 
substrate thicknesses. The losses are particularly great for silicon ground 
planes (0.0025fr-cm) and for small dielectric thicknesses (which result 
in low values of Z ). The losses decrease rapidly as line width increases. 
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The near-end crosstalk for semi-infinite lines is shown in Fig. 20 
as calculated without the influence of the magnetic material. For all 
geometries considered in this paper the near-end crosstalk is below 
4 percent and is typically 1 percent or less. The presence of magnetic 
material would reduce this further. 

For comparison, the transmission parameters of a 50fi microstrip 
transmission line on a conventional 0.625 mm thick Alumina substrate 
are 5,15 given in Table IV. 

3.4 Electrical Results 

Effective transmission line interconnections for chip to chip on a 
layered integrated circuit substrate appear realizable. The major dis- 
advantages are variation of impedance with frequency, appreciable 
delay times (and variation of delay time with frequency causing pulse 
distortion) and significant losses. These problems can be minimized by 



17S 




0.0025 0.0050 0.0075 0.0100 0.0125 0.0150 

DIELECTRIC THICKNESS IN mm 

Fig. 15 — 'Impedance and delay time at 1 GHz for silicon dioxide dielectric on 
silicon ground plane. 
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150 - 



125 - 




N 75 - 



Fig. 16— Impedance and delay time at 1 GHz for 0.0125 mm thick silicon dioxide 
dielectric on silicon ground plane. 



careful design (for example, judicious choices to t m and /*«)• The problems 
of delay times and attenuation (characteristics directly proportional 
to the length of the line) will also be reduced by the much greater chip 
packing densities allowable with substrates of high thermal conduc- 
tivity. Nonetheless, the high losses encountered with Si ground planes 
(along with its lower thermal conductivity) make this material a dubious 
substrate choice* except for possible special cases in which wide lines 
could be used for all portions of the circuitry where attenuation need 
be considered and where the ground returns for high current lines 
(for example, power supply lines) could be deposited as surface metal- 
lization. 

* ThiB configuration (0.0025 fi-cm Si as a ground plane) should not be confused 
with the possible use of high resistivity Si (1500 12-cm) as a dielectric with gold 
surface metallization and a metal ground plane. T. M. Hyltin haj measured the 
dielectric loss in microatrips with high resistivity Si dielectrics. 16 Such a dielectric, 
when used in the configuration discussed in this paper, would contribute an addi- 
tional 0.5 dB/cm attenuation to the attenuation calculated above for aluminum 
ground planes. 
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0.4 



0.2 



0.1 1 

FREQUENCY IN GIGAHERTZ 



Fig. 17 — Impedance and delay time for silicon oxynitride dielectric on aluminum 
ground plane. 

In the case of aluminum ground planes, the losses on narrow (2 mil) 
lines are probably tolerable for present circuitry with the shortened lines 
due to higher packing. For longer lines, power leads, and so on, it may 
be necessary to use wider lines and connect to the beam leads of the 
integrated circuit chips through tapered section "transformers". The 
impedance transition need not be as great as the desired width change 
would indicate if t m could be increased under the wide lines to compensate 
for the decreased L and increased C. 

Many present applications for lower speed integrated circuit logic 
do not require matched transmission lines for all or any of the inter- 
connections. The primary electrical requirement is then the completion 
of interconnects between integrated circuit chips with a minimum 
capacitance to ground. For such applications the narrow conductors 
envisioned in this paper would nearly compensate for the increased 
capacitance per unit length due to reduced t d compared with conven- 
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tional 0.625 mm alumina substrates. Use of SiON dielectric (e fl = 3.4) 
would additionally improve this situation. Resistive signal loss in this 
application should be negligible since typical input impedances are 
1-2K ohms. The possibility of shortened lead length due to higher 
packing densities on high thermal conductivity substrates then suggests 
the possibility of improved electrical performance. 

IV. CONCLUSIONS 

The advantages to be gained from a layered substrate are substantial 
from thermal considerations. The modifications of the electrical prop- 
erties resulting from the layered geometry do not appear to present a 
significant obstacle to the development of such substrates. Indeed 
many improvements in electrical performance appear realizable with 
thin line metallization and transmission line interconnections. Layered 
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Fig. 18 — Impedance and delay time for silicon dioxide dielectric on silicon ground 
plane. 
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Al 2 3 ON Al • ▲ 
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Fig. 19— Attenuation at 100 MHz and 1 GHz. 



substrates with massive metal ground planes will adequately satisfy 
rigid mounting plane requirements, and the physical configuration of 
unbalanced microstrip lines allows the use of reliable, inexpensive 
thermocompression bonding of a beam-leaded integrated circuit chip 
to the substrate. 
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APPENDIX A 



Averaging Technique for Improved Series Conversion 

Theorem: Let T(N) be a Junction 0/ varying sign with zero crossi?igs 

enumerated z„ , 2, /", z, , z,>, /" where i is chosen such that 



dT_ 
dN 



> 0. 



If 



Then 

z. 



E T(N) 

A = Z i 



> 



I 

JV-Z l( - 



£ T(N) 



> 



z, 



E 5TW 



E t(n) + E r(iN0 



E TQi) + E W) 




(32) 



(33) 



(34) 



Fig. 20 — Near end crosstalk between adjacent, parallel L 50fi microstrip lines 
(«- - 0.) 
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Table IV — 50ft Striplines on 0.625 mm Alumina Substrates 



Line width 
Delay time 
Attenuation 
Crosstalk 



0.625 mm 
0.05 naec/cm 
0.048 db/cm 
6.75% 



Prooj: (See Fig. 21) 
Write 



from inequality (33), 



so 



£ T(N) = Si 



£ t(n) > 2 r(jv) 




(35) 

(36) 
(37) 



z-l z L+( z L 




Fig. 21 — Enumeration of extrema and zero crossings. 
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from inequality (32) 

I S, - S i + i | > | S i+2 - S i+a I (38) 

so 

I Si - S i+2 | > | S i+J - S i+3 I; (39) 

similarly 

| 8 M - S i+3 | > | S i+2 - S i+i I (40) 

so 

I Si - S i+2 | > | S i+l - S i+a | > | S i+2 - S i+t \. (41) 

Write 

&, = W. (42) 

Inequality (41) holds for any i chosen in agreement with inequality 
(32) as i -* * . Then 

| Si - W ! > | St+t - W ! > | 8 M - W | (43) 

define 

Ai = Si + 2 8t + l (44) 

but from inequality (43) 

| S, - W | > | S l + 1 - W I 
and from inequality (32) 

S, < S.- +1 
so 

II' - S, > S i+i - W, (45) 

giving 

S, + S I+I < 211' (46) 



so 



OTT' 

Ai < 4" = H\ (47) 



From inequality (44) 

A V + 1 - — I" ■■' ^ MS) 



S.- + ] + <S,- + 2 
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but from inequality (43) 

\s i+1 -w\> \s, +2 - w\ 

and from inequality (32) 

Si+» < o,-+i 
so 

S i+i - W > W - S t+a , (49) 

giving 

S l + 1 + S, +2 > 2W (50) 

so 

9W 

A M >^ = W (51) 

from inequality (47) and (51) 

A i < W < A,- + I (52) 

and the theorem is proved. 

APPENDIX B 

Edge-ground Capacitance for Microstrip 

The capacitance from the edges of the finite thickness strip above 
the infinite ground plane is equivalent to the capacitance problem repre- 
sented in Fig. (22a). This is equivalent to twice the capacitance repre- 
sented in Fig. (22b). The conformal transformation W = sin' 1 (z/t d ) 
results in the configuration of Fig. (22c). The capacitance per unit length 
of this configuration can be calculated by equations (26) and (27) of 
the text. This results in 

+ In (1 + {1 - cxp [-cosh- (1 + /.//„)]}*)]■ 



For t, > t 4 , 



(■<•<}, (l + ^J~ln(2 + j 1 
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(a) 



'J Z- PLANE 



-*| t s U— 2t d — M t 5 k- 



(b) 



V= COSH -1 f 1+3- 



-V2 



Fig. 22 — Conformal mappings. 



(O 



W- PLANE 



V2 



Then 



*-*{.. pl±* 



hi 



1+ : 



/, + 2/.V 



2t d + 2/,/ . J ' 



6o6, Li | 2/ rf + 2/, + f2(/, + /.,)(/, + 2Q]> 



The edge-ground capacitance is then 

2e e R ,_ J2k + 2i. + [2(t d + Q(t rf + 2Q]* \ 



C, = 



In 
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